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ABSTRACT 

We present high-resolution 345 GHz interferometric observations of two extreme lumi- 
nous {L IR k, 10 13 Lq), submillimetre-selected galaxies (SMGs) in the COSMOS field with 
the Submillimeter Array (SMA). Both targets were previously detected as unresolved point- 
sources by the SMA in its compact configuration, also at 345 GHz. These new data, which 
provide a factor of > 3 improvement in resolution, allow us to measure the physical scale 
of the far-infrared in the submillimetre directly. The visibility functions of both targets show 
significant evidence for structure on ~ 0.5 — 1 arcsec scales, which at z > 1.5 translates 
into a physical scale of ~ 5 — 8 kpc. Our results are consistent with the angular and physical 
scales of two comparably luminous objects with high-resolution SMA followup, as well as 
radio continuum and CO sizes. These relatively compact sizes (< 5 — 10 kpc) argue strongly 
for merger-driven starbursts, rather than extended gas-rich disks, as the preferred channel for 
forming SMGs. For the most luminous objects, the derived sizes may also have important 
physical consequences; under a series of simplifying assumptions, we find that these two ob- 
jects in particular are forming stars close to or at the Eddington limit for a starburst. 

Key words: galaxies: high-redshift - galaxies: interactions - galaxies: starburst - infrared: 
galaxies - submillimetre - instrumentation: interferometers 



1 INTRODUCTION 

In recent years, it has become increasingly clear that the 
most infrared(IR)-luminous galaxies pay an important role in 
the star formation history of the Universe. Though they pro- 
vide o nly a trivial contribu ti on in the local Universe, both 
theory ( jHopkins et"al] l2010bl. iHopkins & Hernouistl l2010h and 



Magnelli ct al 



observations i Le Floc'h et al 



2005; Perez-Gonzalez et al 



20091 : iGoto et al.ll2010h have shown that by 



2005; 



z k 1 

the total IR luminosity density is dominated by luminous (10 11 < 
Lir/Lq < 10 12 ; LIRGs), and at higher redshifts by ultralumious 
systems (Lir > 10 12 Lq ULIRGs). Thus detailed study of these 
populations, and the processes driving their tremendous radiative 
output is critical to a thorough understanding of galaxy formation 
more generally. 

Submillimetre-selected galaxies (SMGs), discovered in the 
first deep cosmological surveys at 850/im JSmail et alj 1 1991 



iHughes et alll998l : lBarger et all 19981) by the Submillimetr e Com- 
mon User Bolometer Array (SCUBA: iHolland et alii 19991) . repre- 
sent some of the most extreme obje cts in the high-redshift Universe 
(for a review, see lBlain et al.|[2002h . Owing to a strong negative k- 
correction, the selection function at ~ 800 — 1000/im is flat from 
z ~ 1 — 10, and therefore prov des an unbiased view o f star for- 
mation out to very high redshift dB lain & Longahll 19931) . Though 
their bolometric energy output rivals that of luminous quasars, 
SMGs are several orders of magnitude more numerous at compara- 
ble redshifts {z ~ 2.5: Ichapman et aul2005h . At the same time, 
these two extreme popula tions are thought to be connected via 
an evolution ary sequence jHopkins et al.1 1201 Obi : iNaravanan et al.l 
l20ld.l2009allbl) driving the formation of the most massive galaxies 
jScottetal.N2003 ; iBlain et al.l [20041 : ISwinbank et at] 120061 . 120081 : 
IViero et alj |2009[). Thus, in the context of a merger-driven c os- 
mic cycle (e.g. JSanders et alJI 19881 : IHopkins et alJl2006Ll2008bl Iah. 
SMGs represent the transition objects between star formation and 
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AGN-dominated systems, and as such a compelling laboratory for 
testing models of galaxy formation and evolution in the most ex- 
treme environments. 

While their bolometric luminosi ty is thought to be pr i- 
marily powered by star formation i Alexander et alj 2005bllii , 



6y 

2008l;IValiante et alfcOoilMenendez-Delmestre et alfeOOlbOOS 
Pope et all2008l ; lMomiian et alj201fj| : ISerieant et al.l201Ch . SMGs 
could in principle represent one of two very differe nt channels: a 
stead y -state mode wherein large and extend ed (e.g.. lKaviani et all 
2003; Efstathiou & Rowan-Robinson 2003) reservoirs of gas in 
disk galaxies fuel e d largel y by cosmological g as accretion (e.g., 
iKeres et al.l 12001 l2009bl lal: iDave et all [2009 11 yields a very 
high steady-state star formation rate, or interaction-driven events 
wherein close passages and/or merge rs induce strong bars that 
centrally concentrat e the gas supply dBarnes & Hernquistlll996l : 
Hop kins et ai]|2009h . l eading to a brief and intense burst of nuclear 
star formation (see also Hernquist 1989; Barnes & Hernquist 1991; 
Mihos & Hemquistlll994 Il996l ; fcox et alj|2008t iDi Matteo et"al 



2008b . These two modes can be distinguished by the size of the 



star-forming region, with steady-state disks extending over tens of 
parse cs (R « 28(S„/10 mjy) 1/2 kpc at a = 2: iKaviani et al.l 
l2003h versus interaction-driven bursts which are much more con- 
centrated. There is no a-priori reason to favor either scenario, and 
indeed both very gas-rich d isks (e.g.. lErb et alj|200r3 ; iDaddi et al.l 
2009cj; iTacconi et al.ll2010h and merger-driven, IR -luminous star- 
bursts (e.g jTacconi et alj|2008l ; IWalter et alj|2009l) have been ob- 
served in the high-redshift Universe. 

To date, observational estimates of the physical scale of star 
forming regions in SMGs prefer the interaction-driven scenario. 
However, the majority of detections are in the radio, where by 
leveraging the remarkably tight far-IR radio correlation - in which 
the total rest-frame 20cm luminosity provides a proxy for the to- 
tal star for mation rate among star- forming systems in the local 
Unvierse jCondonlll992l; lYun et alj|200ll; iThompson et alj|2006l; 
lLacki et alj 120091; lLacki & Thompson! l2009h - IChapman et all 
d2004) . lBiggs & IvisorJ d2008h . and lMomiian et alj feoich inferred 
typical phy sical sizes of 5-10 kpc for SMGs. At the same time, CO 



imaging dNeri et alj 12003] ; iGreve et al J 12005); ITacconi et alj l200d 



120081 : Daddi et aTT l2009bll£j: iBothwell et alj boich tends to yield 
comparable sizes. However, each of these techniques has its short- 
comings: the spatially resol ved far-IR/radio corre l ation is still not 
well understood (see e.g., Hipp elein et ail 120031 : iMurphvl 120061 : 
Tabat abaei et alj2 007). and imaging redshifted CO emission in ro- 
tational lines with excitation requirements higher than the J = 
1 — > transition may substantially underestim ate the starburst 
scale dNaravanan et aU2008l : llvison et al]|2009ah . 

Therefore, it is preferable to measure the starburst scale 
in the far-IR directly. However, the beam sizes characteristic of 
single-dish instruments are often nearly an order of magnitude 
too large to probe the relevant spatial scales. And, owing to tech- 
nical constraints, even interferometric imaging is typically too 
coarse. In fact, to date virtually all SMGs with interferometric fol- 
lowup are compact relative to the beam (~ 2 - 3 arcsec FWHM ; 



Dannerbauer et al. 2004, 2008; Iono et al. 2006; Wang et al. 2007; 
Younger et al. | l200l l2008j, l2009bl : ICowie et alj|2009h . Recently. 
Youn ger etal.l d2008b) imaged two of the brightest SMGs known 



1 Here disk refers to a stable, gas-rich disk rather than one that is clump- 
unstable and highly turbulent jElmegreen et al l |2008| ; iDekel et alj |2009| ; 
ICeverino et alj2010h . This scenario is discussed in ij |4.4| in the context of a 
potentially mult-component far-IR morphology. 
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Figure 1. The u — v coverage for our high resolution interferometric imag- 
ing of AzTEC4 and AzTEC8 for both the compact (COM: black), and ex- 
tended (EXT: red) configuration. For further details, including weather con- 
ditions and on-source integration times, see TablefT] 



with the Submillimetre Array (SMA: Ih"o et al ,l2004h in its most ex- 
tended configurations - the highest resolution submillimetre imag- 
ing of high-redshift starbursts achieved to date - and found that 
their far-IR emission was extended on ~few x kpc scales. Such ex- 
tremely bright objects were particularly interesting because at these 
size scales and luminosities the implied energetics start to run up 
against physical limits on t he star formation rate imposed by the 
dynamical time of the gas telmegreerJI 19991) or the effects of ra- 
diation pressure on the dust dThompson et al.ll2005l : iMurrav et al.1 
120051; lMurray||2009l ; iHopkins et all2010al) . 

In this work, we present similar observations of two other 
exceptionally bright, and therefore lumino us objects in a 1 .1mm 
blank-field survey of the COS MOS Field dScott et alj|2008l) per- 
formed with the AzTEC camera dWilson et al]2008l) : AzTEC4 and 
AzTEC8. In both cases, high-resolution follow-up observations 
provide a measurement (rather than an upper limit) on the scale 
of the far-IR emission, and therefore the starburst itself. This pa- 
per is organized as follows: in §[2] we summarize the observations 
and data reduction, in §[3]\ve present the results, in §|4]we discuss 
the implications, and in §[5] we conclude. Throug hout this work, 
we ass ume the 7-year WMAP cosmological model iKomatsu et al.1 
(2010), though any set of cosmological parameters within reason 
will not qualitatively change the interpretation of these results. 



2 OBSERVATIONS AND DATA REDUCTION 

The two targets - AzTEC4 and AzTE C 8 - were selected from the 
SM A/ AzTEC interferometric survey d Younger et al]|200l [2009b) 
of bright source s in the AzTEC 1.1mm (270 GHz) survey of the 
COSMOS field dScott et alj|2008l) . They were first identified in the 
AzTEC map with deboosted flux densities of 5.2^^4 and 5.5^'3 
mjy, and were later detected by the SMA in its compact (COM) 
configuration at 870/im (345 GHz) with flux densities of 14. 4± 1.9 
and 19.7 ±1.8 mjy (assuming a point-source model). We then tar- 
geted both sources with the SMA in its extended (EXT) config- 
uration with a goal of measuring the source size. The observing 
conditions for and an overview of all available data for the two tar- 
gets are summarized in Table [T] and the u — v coverage presented 
in FigureQ] 

A detailed description of the calib ration strategy for the COM 
configuration tracks are provided in lYounger et alj d2007l) and 
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Table 1. Track Details 



Target 


Configuration™ 


u — v Coverage 


Beam Size 


Date 


< T225GHz > 


Obs. Time 6 


Reference 






[kA] 


[arcsec] 


[dd.mm.yy] 




[his] 




AzTEC4 


COM 


20-75 


2.71 X 2.12 


21.01.07 


0.06 


5.9 


Y07 




EXT 


50-250 


0.86 x 0.77 


23.02.09 


0.05 


3.9 


This work 


AzTEC8 


COM 


20-75 


2.69 x 2.19 


17.12.07 


0.05 


6.2 


Y09 




EXT 


50-250 


0.86 x 0.55 


16.02.09,22.02.09 


0.1,0.06 


3.7,3.8 


This work 



a For details on SMA configurations, see |http : / / smal ■ sma ■ hawaii ■ edu/status ■ html ° T otal on-source integration time in that configuration. ' 

Y07: lYoungeretal.l J2007t) ; Y09: Youn ger etai J <2009bl) 



870/im 



Band 



Ks 



20cm 



AzTEC4 




AzTEC8 








Figure 2. Postage stamps for AzTEC4 (top) and AzTEC8 (bottom) including (left to right): 870/^m SMA (natural weighting), Subaru B-band, HST/ACS 
i-band, CFHT Ks-band, and VLA 20cm imaging. For both sources, we only include SMA data from the EXT tracks. The red contours indicate 3, 5, 7, ... X 
the r.m.s. noise level, the grey hashed ellipse the beam size, and the grey dashed line 1 arcsec for scale. 



lYounger et a l. (2009b) for AzTEC4 and AzTEC8 respectively. For 
the EXT tracks presented in this work, the receiver was tuned to 
340 GHz in the LSB, and averaged with the USB for an effective 
bandwidth of 4 GHz cen tred at 345 GH z. Passband calibration was 
performed using 3C84 teennetd[l962h . and primary flux calibra- 
tion was done using either Ceres (on 16 and 23 Feb 2009) or Titan 
(on 22 Feb 2009). The target was observed on a 10 minute cycle 
with two primary gain calibrators: 1058+015 (~ 2 Jy; 15 degrees 
away) and 0854+201 (~ 6 Jy; 24 degrees away). Because Ceres 
is known to be variable at the ~ 20 — 30% level due to rotation 
dAltenhoff etalj 1 1994 iRedman et al] 1 19981: iBarrera-Pineda et all 
12005b . we confirm this flux scale by checking that the flux density 
for 0854+201 derived using Ceres as the primary flux calibrator 
(6.0 ± 0.3 and 5.7 ± 0.3 Jy) were consistent with those using Titan 
(5.7 ±0.3 Jy). 

In addition to the two primary targets, we observed a nearby 
test quasar - J 1008+063 (~ 0.2 Jy) which was 5 degrees away 
from the targets - once every 60 minutes throughout the track to 
empirically verify the phase transfer and inferred source structure, 
and estimate the systemati c positional uncertainty. This source is 
included in both the JVAS (Patna ik et al] 19921; [Browne et aljl998b 
and VLBA Calibrator dMa et al.ll998l;lBeaslev et al.ll2002l) surveys 
of compact, flat-spectrum radio sources, has an absolute position 
known to be tter than 20 mas, and w as confirmed to be compact at 
345 GHz by I Younger etai] d2008bb . 



We also make u se of extensive mult iwavelength data in the 
COSMOS Field (see IScoville etalj l200l for an overviewf] in- 
cludin g: Subaru ground bas ed optical and near-IR l Taniguchi et al.l 
|2007|) . HST/ACS i-band (Koekemoer etai] |2007|). IRAC 3.6- 



P 



8/xm and MIPS 24u m dSanders et al] |2007|). VLA 20cm 
( Schin nerer et al.l 120071) . and Chandra X-ray dElvis et al] 12009]; 



Puccetti et al.ll2009l) imaging. 



3 RESULTS 

Imaging and model fitting for the calibrated v isibility data wa s 
performed using the MIRIAD software package dSault et ai]|l995l) . 
We utilize a natural weighting scheme, yielding dirty maps for the 
EXT tracks that had r.m.s. noise levels of 1.6 and 1.8 mjy/beam 
for AzTEC4 and AzTEC8 respectively. The peak flux density in 
each map was spatially consistent with detections in earlier COM 
tracks and highly statistically significant: for AzTEC4 the peak was 
8.0 mjy/beam (S/N-5), and for AzTEC8 it was 12.8 mjy/beam 
(S/N-7). 

Postage stamps of these maps, along with contours overlaid 



2 More current information is 

http : //cosmos . astro. caltech.edu/ 



available 
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Figure 3. Binned and vector-averaged real visibility amplitudes as a function of u — v distance for the target (top row) and test quasar (bottom row). The 
targets - AzTEC4 (left) and AzTEC8 (right) - and the calibration strategy are discussed in §|2] For the targets we include all available data, including both 
COM and EXT tracks, details of which are provided in Table[T]and the u — V coverage presented in Figure[T] for the test quasars we only include the EXT 
data. For comparison, we provide a series of symmetric Gaussian source models with FWHM sizes of (left to right) 10, 5, 2, 1, 0.5 arcsec. The flat solid and 
dotted lines indicate a point-source fit (including only the COM data for the targets, and only the EXT data for the test quasars) and the associated statistical 
uncertainty. We find that both AzTEC4 and AzTEC8 show clear evidence of structure on 0.5 — 1 arcsec scales, and the test quasars confirm the phase 
transfer. 



Table 2. Positions and Source Structure 



Name 


Config." 


Model 


a 


<5 




A<5 6 




e c . 

maj 


e c . 

mm 










[J2000] 


[J2000] 


[arcsec] 


[arcsec] 


[mJy] 


[arcsec] 


[arcsec] 


[deg] 


AzTEC4 


C 


Point 


09:59:31.72 


+02:30:44.0 


0.15 


0.24 


14.4 ± 1.9 










C+E 


Gaussian 


09:59:31.709 


+02:30:44.06 


0.11 


0.08 


13.1 ± 1.8 


0.6 ±0.2 


0.4 ±0.2 


Q 




C+E 


Disk 


09:59:31.709 


+02:30:44.06 


0.09 


0.07 


13.1 ± 1.7 


1.0 ±0.4 


0.7 ±0.6 


30 


AzTEC8 


C 


Point 


09:59:59.34 


+02:34:41.0 


0.10 


0.10 


21.6 ± 2.3 










C+E 


Gaussian 


09:59:59.334 


+02:34:41.12 


0.064 


0.058 


17.7 ±2.3 


0.6 ±0.2 


0.5 ±0.3 


40 




C+E 


Disk 


09:59:59.334 


+02:34:41.09 


0.054 


0.056 


17.2 ± 2.1 


1.0 ±0.5 


0.4 ±0.8 


20 



a See Table[T]and Figure[T]for details. b Combined statistical and systematic uncertainty, where the systematic uncertainty is estimated from the position of 
the test quasar. c f m aj and $ m j n represent the FWHM or diameter of the major and minor axes for the Gaussian and elliptical disk models respectively. d 

Position angle. 



on optical, near-IR, and radio data are shown in Figure E] Nei- 
ther source has a statistically significant detection in the optical 
or near-IR, though there is a marginal ~ 23.3 mag K-band de- 
tection associated with AzTEC4. In the radio , as expected there is 
no emission coincident with AzTEC4 (as in lYounger et alj|2007h 
and the centroid of AzTEC8 is consistent with that of an 89 ± 11 
[iJy s ource - as in the CO M data (which has a comparable noise 
leveh lYounger et alj200 9b) the nearby secondary radio counterpart 
of AzTEC8 is not detected in the EXT image. 



We then combined the new EXT with existing COM data and 
inspected the visibility function for each target (see the top panel of 
Figure[3}. When binned by u— v distance and vector-averaged, both 
AzTEC4 and AzTEC8 show significance evidence for structure at 
> 100 kA, which corresponds to angular scales of < 0.6 arcsec. 
This is consistent with the images, which exhibits lower peaks in 
the higher resolution maps: for AzTEC4 the peak flux density in 
COM 13.2 ± 1.7 mjy/beam versus 8.0 ± 1.8 mjy/beam in EXT, 
and for AzTEC8 the peaks are 19.7 ± 1.8 mly/beam in COM and 
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12.8 ± 1.8 mjy/beam in EXT. In both cases, the test quasar em- 
pirically confirms that this structure is real, and not the result of 
decorrelation on the longest baselines (see bottom panels of Fig- 
ure O. 

We then fit the calibrated visibility data to one of two simple 
models - a 2-dimensional Gaussian or an elliptical disk - the results 
of which are presented in Table|2] The best-fit parameters for both 
sources show structure on ~ 0.5 — 1 arcsec scales at > 2<r confi- 
dence. We also confirm that this size measurement is not artificially 
imposed by differing flux scales between the COM and EXT data 
by fitting the same source model to each set of calibrated visibilities 
independently: when we assume e.g. an elliptical Gaussian we find 
total flux densities (for AzTEC4 and AzTEC8 respectively, and in- 
cluding an estimate of the systematic uncertainty in the flux scaljfb 
of 14.1 ± 2.5 sta t ± 0.8 sys and 17.9 ± 3.7 stat ± 0.9 sys mJy for the 
COM data, and 9.2 ± 3.1 sta t ± 0.5 sys and 21.0 ± 4.8 8tat ± 1.0 sys 
mJy for the EXT data; that these are consistent between configura- 
tions to within the statistical and systematic uncertainties indicates 
that the size measurements listed in Table[2]are robust. 



4 DISCUSSION 

4.1 The physical scale of the far-infrared 

These angular scales are qualitatively consistent with size mea- 
sure ments derived from previous high-resolution far-IR imag- 
ing dYounger et al.|[2008bl). as well as other wavelengths includ- 
ing radio continuum dChapman et alj|2004l; iBiggs & rvison||2008| ; 
Momi ian et alj201Ct) and CO dNeri et al.l2003l ; lTacconi et alj2O06L 



2008) imaging. In our view, imaging in the far-IR directly is the 



most robust of these methods because it directly probes the ob- 
scured starburst; though the galaxy-wide far-IR/radio correlation is 
thought to app ly - at least in an qualitative sense - at redshifts typ- 
ical of SMGs dGarretfc002l;lGnippioni et alJl2003l:lAppleton et all 



20041: iBovle et alJl2007lTYounger et al.1 2009a ; 



Ivison et al 



Murphvl 120091) 7 the spat all y resolved correlat ion 



2009b; 



is not 



par- 



ticularly well underst ood dHippelein et al.1 120031 ; iMurphvl 120061 ; 
iTabatabaei et~afll2007t) . while CO imaging is largely restricted to 
excited states with higher critical densities than the bulk of the star- 
forming gas, and there f ore may be biased towards smaller sizes 
dNaravanan et al.l [20081 ; Ifvison et alj [2009a). However, given the 
agreement between these three methods - now including twice as 
many objects with sizes measured in the far-IR directly - we have 
newfound confidence that the general conclusion that the most lu- 
minous SMGs have typical angular scales of ~ 0.5 — 1 arcsec is 
robust. 

To translate this into a physical size requires knowledge of the 
angular diameter distance. Fortunately, at 1.5 ^ z 5 and assum- 
ing th e 7-year WMAP ACDM cosmological model ( Komats u et al.1 
|2010|) . this quantity has a very weak scaling with redshift. In gen- 
eral, the physical size (£) is approximately: 



11 



V arcsec / 



z) ' 20 kpc. 



(1) 



Therefore, existing size measurements for SMGs suggest a typical 
physical scale of ~ 4 — 8 kpc at z ~ 2 — 5. W hile there is clear ev - 
idence for > 20 kpc sizes among some SMGs dlvison et alj2 009a). 



3 Estimated from measurements of the flux density of 0854+201 
which are available online as part of the SMA Calibrator List at 

|http : / / smal . sma . hawaii . edu/ callist/callist. html| 



the data presented in this work strongly argue against pure steady- 
state s tar form a tion fed by cosmo logical gas accretion flKeres et al.l 
I2005L I2009bl l3; lDaveetal.l[2009Ti as the primary mechanism for 
producing SMGs, and rather favors the merger-driven scenario 
dNaravanan etai]|2010l . |2009a) in which interaction-driven star- 
bursts give rise to these extreme objects. This in turn fits into a 
more general evolutionary p icture in which SMGs are transition ob- 
jects (Hopkins et al. 2010b) between isolated, gas-rich spirals and 
lumi nous quasars, and the progenitors of passiv e, elliptical galaxies 
(e.g.. ISanders et al.lll988l ; Hopkins et alj|200rj , l2008tal fah. Further- 
more, since thes e sizes are considerably l arger than local ULIRGs 
(see review by I Sanders & Mirabell 1 1996T) . which themselves are 
powered primarily by major merge rs (e.g., IScoville et "al] |2000l : 
lKimetal.l (2002; Veil leux et alj |2002). they are consistent with the 
canonical picture of SMGs as 'scaled up' ULIRGs. 



4.2 Are the most luminous submillimetre galaxies 
Eddington-limited starbursts? 

Size measurements are particularly interesting for the most lumi- 
nous SMGs because their star formation is sufficiently rapid to 
run up against fundamental physical limitat ions. Owing to a strong 
negative fc-correc tion in the submillimetre dBlain & Longaigl 1993c 
lBlainetal]|2002l) . we can infer the total infrared luminosity in 
an approximate sense by assuming an Arp 220 template, which 
for z > 1 yields Lj/fl/lO 12 L s » 2 x [S 8 70 M m/mJy] (see 
a lso iNeri et alj|2003l ; [Younger et alj |2008a). By further assuming 
a lSalpetej dl955t) IMrrT this can be converted into a star formation 
rate SFR « 340 x [<S l 87u um/mJyl Mq yr" 1 . At the same time, 
Eddington-type arguments (Murrav et al. 2005) - under a series of 
simplifying assumptions such as a spherically symmetric geome- 
try, an isothermal sphere density str ucture, a small volume filling 
factor for molecular gas, and again a Salpete^ dl955h IMF - yield a 
maximum star formation rate of: 



SFR„ 



9000-looAtpcfnoo Mq yr 1 



(2) 



where a is the line-of-sight gas velocity dispersion in units of 400 
km s -1 , K100 is the opacity in units of cm 2 g -1 (usually taken to 
be w l; lMurrav et al.l2005l ; Tfhompson et al.l2005t) . and D kpc is the 
characteristic physical scale of the starburst - either the Gaussian 
FW HM or the elliptical dis k diameter. 

lYounger et ail d2008bl) found that GN20 and AzTEC 1 - two of 
the most luminous SMGs known - were potentially at or close to 
this Eddington limit. In this work, we measured sizes for two com- 
parably luminous objects: taking the flux densities inferred from 
the EXT+COM data, the SFR « 4400 Mq yr" 1 and « 6000 
Mq yr -1 for AzTEC4 and AzTEC8 respectively. Both have char- 
acteristic sizes of I ~ 5 — 8 kpc - where the range represents 
uncertainty arising from a Gaussian versus disk source structure - 
yielding a maximal star formation rate of SFR max ~ 3600 — 7200 
Mq yr -1 . Therefore, both objects appear to be forming stars at or 
near the Eddington limit for a starburst. 



4 As noted in Younger et al. 1 2008b), The inferred luminosity is known to 
be uncertain b y a factor of ~ 2 — 3 due to variations in the dust temperature 
and emissivity Blain et al. 12003). For example, if we assume a Mrk 231 
template with significant AGN contribution to the farlR, the i nferred farlR 
luminosity and SFR will be 



Stevens et al. 2005 



40 - 60 % lower (e.g. 
iHuang et alj2 007). Furthermore, using a lKroupa) d200ll) o rlChabrieJ d2003 ) 
IMF will tend t o lower any SFR by m 40% <Kennicuttlll998UBelllhoOS 
iBell et alj2005l) . 
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There are, however, a number of important caveats to con- 
sider: 

1. Uncertainties in estimating the star formation rate: There 
are few observational constraints on the stellar IMF at high- 



Table 3. Positions and Source Structure 



redshift and in extreme environments (e.g.. | Fardal et alJ 2007: Dave 
2008l:lvan Dokkurrj2008l ; lBaugh et al]2005r . lSwinbank et al.l2008l ; 
Tacconi et alj 120081) and t he shape of the far-IR SEP as a func- 



tion of luminosity ( e.g.. bale et alj l200ll : bale & Helot] 120021 : 
ICharv & Elbaj|200ll ; lLagache et alj|2003ri . 

2. The volume filling factor of molecular gas: In the optically 
thick limit - in which the molecular gas has a volume filling 
factor near unity - the Eddington limit is an order of magni- 
tude higher. A number of rece nt studies have found that this limit 
descr i bes even local ULIRGs dScoville et al.lll99ll: iDownes et al.1 
Il993l : iDownes & Solomon] 1 19981 : ISolomon et alj|l997h and there^ 
fore we might expect it to apply well to the much more extreme 
environments in SMGs. This does, however, require a rather low ef- 
fective dust temperature; in the optically thick limit, the brightness 
temperature at 345 GHz (Tb) is equal to the physical temperature 
(T d ). While the sizes and flux densities of AzTEC4 and AzTEC8 
require Td ~ 5(1 + z) K (see Figure 5 and discussion in § 4 of 
lYounger et al1 l2008b). which at z ~ 2 — 4 is rather low compared 
to other SMGs for which Td has been measured independently via 
far-IR SED fitting « T d >~ 35 K for typical SMGs, and the most 
luminous objects m ay be preferentially hotter: iKovacs et al ] |2006l : 
Coppin et al]l2008h . it is certainly not out of the question. 

3. The depth of the potential well: The maximal star formation 
rate has a strong dependence on the gravitational potential, which 
enters as erJoo • Dynamical masses measured via CO s pectroscopy 
have found that, in the m ean 0400 ~ 1 for SMGs I IGreve et al.1 
120051 : ITacconi et al] [2006). However, merger-driven models for 
SMGs sugg est that the brightest 850/xm sources are also the most 
massive dNaravanan et alj|2010h . a nd therefore would likely have 
0-400 Z 1.5 (for GN20, 0-400 ~ 1.4; ICarilli et alj2010h or so which 
could lead to significantly higher Eddington limits, even in the op- 
tically thin limit. 

Generally speaking, however, we find that AzTEC4 and AzTEC8 
are extende d on physical scales com parable to the other two SMGs 
presented in l Younger et alJd2008bl) . and are potentially radiating at 
or close to the Eddington limit for a starburst. 



4.3 Obscured AGN versus star formation 

Though in gen eral SMGs are thoug ht to b e star forma- 
tion d ominated ([Alexander et all l2005blaL 120081: IValiante et alj 
2007; Menendez-Delmestre et al. 120071 120091 ; IPope et alJ 12008; 
Momiian et alj|20ld : ISerieant et alJl2010I) . in principle they could 
contain a significant contribution from an obscured AGN. In fact, 
a preliminary analysis of recent very long baseline interferometry 
observations provides strong evidence that the radio continuum in 
at lea st some SMGs is powered primar ily by an ultracompact AGN 
core dBiggs et alj2009l ; lYou nger 2010). Clearly if this was the case 
for either AzTEC4 or AzTEC8 it would severely compromise their 
interpretation as extreme starbursts at or near their Eddington limit 
- the Eddington limit for a ~ 10 9 Mq supermassive black is well 
in excess of 10 13 Lq. 

Re cent X-ray imagin g of the COSMOS field (the C-COSMOS 
Survey: lElvis et alj|2009l) provides some constraints on the AGN 
content of these objects; though starbursts also produce significant 
X-ray emission, a detection in the hard band (2 — 8 keV) would 
be strong evidence for the presence of a buried AGN - particularly 



Name 


pa 
1,890/jm 


pb 
2,S90Mm 


A6» c 




[mJy] 


[mJy] 


[arcsec] 


AzTEC4 


5.8 ±2.5 


6.8 ±2.6 


0.48 ±0.12 


AzTEC8 


13.0 ±1.6 


5.4 ± 1.6 


0.75 ± 0.19 



a Flux of the first component derived from fitting a two-component point- 
source model to the calibration visibilities. b Flux of the second component 
derived from fitting a two-component point-source model to the calibration 
visibilities. c The separation of the two fitted components. 

at high-redshift. Therefore, we have examined the C-COSMOS X- 
ray imaging data for AzTEC4 and AzTEC8. While AzTEC8 shows 
no evidence for a detection, AzTEC4 exhibit s a tentative hard X- 
ray source. A formal source extraction (as in IPuccetti et aUl2009l) 
yields net counts in the hard band of 5.5 ± 2.7 cts, which translates 
into a flux of Fhx = 8.0 ± 3.9 erg cm -2 s -1 - a ~ 2a detection . 
If we assume the bolometric corrections of Hop kins et al. I d2007h , 
this implies an AGN with bolometric luminosity L oa i /10 11 Lq = 
1.2 ± 0.5, 3.4 ± 1.4, and 7 ± 2.4 at z = 2, 3, and 4 respectively, 
which translates into M B Hrj a dd/W 7 Mq = 0.4 ± 0.2, 1.0 ± 0.4, 
and 2.1 ± 0.8 where r/ e dd is the Eddington ratio of the AGN (typ- 
ically near unity during the peak of the starburst: Hopki ns"et al.1 
l2005tal R>. While we cannot rule out significant X-ray absorption 
(Compton-thick in the case of non-detections), at these redshifts 
and energies the optical depths are unlikely to be extreme. There- 
fore, the X-ray data suggests that AGN do not contribute signifi- 
cantly to the IR luminosity of AzTEC4 or AzTEC8. 



4.4 A multi-component source structure? 

The visibility functions for AzTEC4 and AzTEC8 show clear ev- 
idence for structure on ~ 0.5 — 1 arcsec scales (see Figure [3}, 
and fitting a model to the visibility data yields a statistically sig- 
nificant size measurement (see Table [2}. However, these observa- 
tions are not sufficiently high-resolution, nor do they have suffi- 
cient signal-to-noise to distinguish an extended source structure 
from multiple compact components separated by less than the beam 
size. When we fit a dual point-source model to the visibility data, 
both sources - especially AzTEC8 - yield statistically significant 
(> 2a) measurements for the implied sub-component flux densi- 
ties and separations (see Table [3} . Therefore, while the visibility 
functions are consistent with an extended source structure, the data 
could also indicate a multi-component source structure produced 
by either more compact starbursts or Compton-thick AGN. How- 
ever, the sizes listed in Table [2] can be thought of as upper limits 
on the physical scale of the starburst in these objects, and therefore 
still argue against an ext ended star-forming disk. 

ICarilli etai] d2010h h ave argued - using the spec ific case of 
GN20 (see also lPope et aUl200d : lYo~unger et alj|2008bh - that such 
a multi-component source structure, particular one that appears 
circularly supported with significant extra-nuclear star formation 
is inconsistent with a major merger; rather, they they find that 
it is indicative of clumpy star format ion owing to local in stabili- 
ties in turbulent, high-redshift disks dElmegreen et alj |2008). This 
is, however, not the c ase. First, both hydrodyna mical simulations 
dCeverinoet"aH f2010) and analytic arguments dDekel et"al] [2009) 
indicate that this mode of star formation is steady-state, while 
the specific star f ormation rate (SSF R) of GN20, for example, is 
» 10-30 Gyr^ 1 dCarilli et al.l2010h : a long duty cycle at this high 
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a SSFR is unphysical, and furthermore is inconsisten t with cold- 
mode accretion rates from cosmological simulations I Ker es et"al] 
l2009bl fah. Second, significant non-nuclear CO emission is not in- 
consistent with the merger scenario, particular one involving more 
than two participatin g galaxies (though ad mittedly a more ex- 
treme case, see e.g. iNaravanan et alj|2006l) . Finally, should the 
clumps appear to exhibit disky kinematics, this would also be con- 
sistent with such a gas-rich merger; hydrodynamical simulations 
have shown that disky gas kinematics are generically conserved 
when the progenitor s are extremely gas-rich (Ro bertson et alj2004l ; 
iHopkins et aU2009l) . 



5 CONCLUSION 

We present high-resolution interferometric observations of two 
SMGs in the COSMOS field - AzTEC4 and AzTEC8 - performed 
with the SMA at 345 GHz. The targets, two of the most luminous 
SMGs known, were previously detected as compact sources with 
the SMA at the same fre quency in its COM configu ration (beam 
size « 2.7 x 2.2 arcsec: I Younger et al .1 [20071. 1 2009bh . These new 
observations, which offer a factor of ^ 3 improvement in resolu- 
tion, allow us to measure the angular size of the two targets. The 
visibility functions show significant evidence of structure on angu- 
lar scales of « 0.5 — 1 arcsec, in agreeme nt with the sizes of tw o 
comparable objects measure d by the SMA Younger et alj[2 008b), 

~~ ™ J/-*! ~ — „«. „i or\n/il.lr>: t.,: 1 



as well as radio continuum ^Chapman et al. 2004 ; Bi ggs & IvisorJ 
2008| : lMomiian et alj201oh and CO jNeri et al.l2003l ; lTacconi etaU 
20061 120081 : iBothwell et ai]|2010h imaging. Owing to the weak 
scaling of angular diameter distance with redshift for 2 1.5, we 
can convert this angular scale to physical units and find that the 
far-IR in these SMGs is extended over ~fewx kpc. This provides 
evidence in favor of a merge r-driven scenario for forming SMGs 
1 Narayanan et al. 2009a, 2010), rather than extended gas-rich disks 



Dave et al.ll2009l : llvison et al. 2009a). For the most luminous ob- 



jects, the derived sizes may also have important p hysical conse- 
quences ; as with two comparable objects studied by Young er"et al.l 
(2008b), under a series of simplifying assumptions the two targets 
in this stud y are forming stars at or near the Edding ton limit for a 
starburst dMurrav et alj2005l ; lThompson et alfeooa) . 
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